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AbSTRACT 


An analysis and a numerical lifting surface method are 
developed for predicting the unsteady airloads on two-dimen- 
sional circulation control airfoils in incompressible flow. 

The analysis and the computer program are validated by corre- 
lating the computed unsteady airloads with test data and also 
with other theoretical solutions. 

Additionally, a mathematical model for predicting the 
bending-torsion flutter of a two-dimensional airfoil (a reference 
section of a wing or rotor blade) and a computer program using 
an iterative scheme are developed. The flutter program has a 
provision for using the CC airfoil airloads program or the 
Theodorsen hard flap solution to compute the unsteady lift and 
moment used in the flutter equations. The adopted mathematical 
model and the iterative scheme are used to perform a flutter 
analysis of a typical CC rotor blade reference section. The 
program seems to work well within the basic assumption of the 
incompressible flow. 



NOTATION 


C. jet blowing momentum coefficient 

D semi-width of a panel 

h,H position of the wake jet 

k,K doublet strength 

i semi-chord length of an airfoil 

p frequency of oscillation 

t , T time 

x,X coordinate along the chord line of an airfoil 

U free stream velocity 

w,W downward displacement of an airfoil 

Z coordinate for plunging motion 

a angle of attack, coordinate for pitching motion 

Y,r vorticity 

a) reduced frequency 

P density 
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INTRODUCTION 


In the last few years there has been a considerable in- 
crease in activity in the prediction of unsteady airloads on 
the circulation-control and jet flap airfoils. Williams et 
al. (Ref. 1) and Platzer (Ref. 2) presented a review of the 
experimental and analytical studies of jet flap airfoils. 

Several investigators conducted a series of experimental studies 
on airfoils with trailing-edge jet flaps. Simmons and Platzer 
(Ref. 3) measured the frequency response of a trailing-edge 
flap on a two-dimensional NACA 0012 airfoil. Two values of 
jet momentum coefficients (Cj = 0.14 and 0tl58) and several 
values of reduced frequencies were investigated at constant jet 
deflection amplitudes and the airfoil at zero angle of attack. 
Trenka (Ref. 4) considered the jet flap and wing oscillations 
on a two-dimensional NACA 0015 airfoil with a 20% flap over the 
span. Takeuchi (Ref. 5) measured the oscillatory forces and 
moments on a semi-span wing model of a NACA 632A015 airfoil 
with a trailing-edge jet flap. Using false walls, various 
aspect ratio wings ranging from 2.1 to infinity were simulated. 
Kretz (Ref. G) investigated the harmonic oscillations of a jet 
flap on a NACA 0012 at Cj » 0.25, 0.50 and 1.0 for a reduced 
frequency range from 0 to 1.26. 

Some progress has been made towards the development of 
theoretical solutions for jet flaps under unsteady flow condi- 
tions using potential flow theory. It has been a standard pro- 
cedure to assume pure harmonic motions of jet and airfoil. 

Three types of motions are of interest; namely: (i) the sta- 

tionary airfoil with oscillating jet flap; (ii) the pitching 
motion of the airfoil; and (iii) the plunging motion of the 
airfoil. For harmonic motions of the airfoil or jet flap, the 
pressure distribution, and hence the aerodynamic loads, lag the 
motion. The prediction of this lagged response taking the jet 
interaction into account has not been too successful even under 
potential flow conditions. Spence (Ref. 7) suggested that for 
low values of reduced frequency (less than 0.2), the wall jet 
effects may be neglected and Theodorsen's trailing-edge, hard- 
flap solution may be a good blowing representation. At low 
blowing coefficients, there is a reasonable correlation between 
the experimental and analytical (Theodorsen) results for isolated 
airfoils. Theodorsen's solution can easily be modified to in- 
clude the wake effects of a helicopter rotor blade in hover using 
Loewy ' s mathematical representation of a wake model. 
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Spence* (Ref 7) obtained a solution for jet flaps in un- 
steady flow, modifying the governing flow equation and the jet 
dynamic equation using a weak jet approximation. However, his 
solution is valid only for small values of jet momentum coeffi- 
cient. Potter (Ref. 8) solved the problem as an initial value 
problem by reducing the governing equation into a finite-dif- 
ference form and using point vortex distributions. The jet 
is initially specified as a straight horizontal line with zero 
vorticity, after which time the jet exit angle is either in- 
stantaneously deflected or continuously oscillated, and the 
shed vortices are assumed to propagate downstream with the 
f ree-st ream velocity. Conservation of circulation is maintained 
tor each time stop, and the jet shape is determined by taking 
into account .ill changes in bound airfoil vorticity in comput- 
ing the downwash and the resultant velocities at each jet vor- 
tex point. The adoption of an averaging scheme over several 
neighboring vortices provided a satisfactory approach to cor- 
rect the erratic vortex motion with increasing time steps, and 
gave satisfactory results for jet momentum coefficients smaller 
t han 0.1. 

Scliniidt (Ref. 9) presented a complete summary »nd evaluation 
ot t lie experimental investigations (Refs. 10 through 13) on 
the circulation control airfoils under unsteady flow conditions 
conducted at the Naval Postgraduate School. Bauman (Ref. 10) 
developed and tested the operation of a control valve in the 
pressure supply lino just upstream of the airfoil cavity so as 
to permit the superposition of a harmonic perturbation upon a 
mean value of cavity pressure. Rickelsimer (Ref. 11) developed 
a numerical algorithm for ident ifying Fourier components in 
multi-channel, discretized, truncated signal information. This 
technique enabled the determination of amplitude and phase 
ngle of the various data channels relative to the reference 
clock channel (unsteady cavity pressure) and also the in-phase 
and out-of-phase lift coefficient components from unsteady pres- 
sure data. Kail (Ref. 12) investigated the behavior of the 
Coanda sheet formed by the tanqential jet flow over the airfoil's 
rounded trailing edge and noted that the superposition of an 
oscillatory blowing component upon the steady tangential jet 
flow did not produce any discernible change ir the average 
value of lift augmentation. Englehardt (Ref. 13) developed the 
hardware and software implementation of a microprocessor based 
on a high-speed digital data acquisition and reduction system 
that was tailored for use in time-varying signal analysis of 
he unsteady pressure information derived by oscillating the 
cavity pressure of a circulation-control airfoil. Schmidt 
l Ret 9) analysed the results for the CCA model at one angle 
of attack, one tunnel airspeed, with the harmonic oscillation 
of the cavity pressure at four frequency values using one 
average value of momentum blowing coefficient, and arrived at 
the following conclusions: 
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a. The airfoil behaved in a linear manner. 

b. The mean or average value of the lift coefficient was not 
altered during unsteady blowing. 

c. The frequency response traits of the Coanda sheet region 
suggest the presence of a transportation lag between the 
pressure cavity and the blowing slot. 

d. The overall airfoil lift transfer function from harmonic 
circulation-control variations had a behavior quite similar 
to that of a simple pole in classical control theory. 

e. A frequency-dependent pitching moment about the 54 percent 
chord center of pressure location for circulation control 
type lift was evident with the attributes of aerodynamic 
damping being shown. 


Williams et al (Ref. 1) presented a discussion on the 
appropriate mathematical modeling of the circulation-control 
airfoil under unsteady flow conditions. With a rate of change 
of angle of attack, heave or velocity at constant duct pressure, 
they have noted that the airfoil will produce a corresponding 
lagged response in pressure distribution and lift and have 
hypothesized that these unsteady effects may be treated by 
neglecting the wall jet dynamics, especially for low values 
of reduced frequency. Under this assumption, the CC airfoil 
problem can then be stated in terms of a bluff, elliptical type 
of airfoil with a small trailing-edge vane or flap, a variation 
of whose position produces a corresponding change in circulation 
on the airfoil and vorticity in the wake. For oscillating 
angles of attack or heave, the flap position is assumed to be 
controlled cnly by the blowing coefficient so that the problem 
can be treated classically. Also presented in Reference 14 
was some limited test data which demonstrated that the oscil- 
lating mechanical flap is a good blowing representation. 

In conclusion, at the present time, there do not exist any 
analytical and/or computational methods for predicting the 
unsteady airloads on circulation-controlled airfoils for high 
blowing coefficients. For the present development, a numerical 
lifting surface-method will be adopted. Rao and Jones have 
successfully applied this method for predicting steady and un- 
steady airloads over a variety of thin lifting surfaces in sub- 
sonic flow, such as helicopter blades, cascade of blades, and 
delta wings (Refs. 14, 15, 16 and 17). In all of their work, 
they have used potential flow and small perturbation assumptions 
at small incidence angles. For circulation-control airfoils, the 
effect of the unsteady jet interaction can be incorporated by 
prescribing the jet displacement taking into account the jet 
dynamics . 
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in the present study, a numerical lifting-surface method 
is developed for predicting the unsteady airloads on circula- 
tion-control airfoils, for arbitrary combinations of pitching 
and plunging motions, longitudinal velocity variations, and 
variations in blowing pressure and jet deflection angle. The 
analysis and the computer program are tested by correlating the 
computed unsteady airloads with test data and the approximate 
prediction based on Theodorsen' s method. 

In addition, a simple two-dimensional bending-torsion 
flutter analysis and a computer program are developed using 
the unsteady airloads program developed in the present study. 
This analysis and the computer programs are used to perform a 
flutter analysis of a typical CC rotor blade section. 


UNSTEADY AIRLOADS ANALYSIS 


T he Governing Flow Equations 

The governing equation for incompressible flow relating the 
dow.nwash on the airfoil to the vorticity distribution from the 
leading edge to infinity and the jet dynamic equation relating 
the jet displacement to the local vorticity distribution (see 
Figure 1) are expressed in a ncn-dimensional form as: 


W ( X , T) 



( 1 ) 


and 


(iw + 


±) 

dX ; 


“ - 2y 


d 3 a. 
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( 2 ) 


where 


x = x/i, t = u t/t, w = w/u, r = y/v, h ■ h/i' “* ° p£/u. 

p = r/4, and i is the semichord. 
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For an oscillating airfoil or jet it is assumed hj » 0, and 
h'j * h'j c * s constant at the trailing edge of the airfoil. 

These boundary conditions can be matched by amuiming t lie following 
form of the wake shape: 
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Numerical Lifting-Surface Method 


The numerical procedure uses doublet distributions on the 
airfoil surface. Hence Eqns. (1) and (2) are rewritten in terms 
of doublet distributions as: 


W (X, T) = 
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where 



or 


r (X) = dK/dX, 


and K. (Kj(X)ex =1) is the doublet strength at the 
trailing edge of the airfoil. 


Equations (3) and (5) are combined to yield 
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The last expression in Eqn. (7) can be analytically integrated 
for any given value of X. The expansion of this wake integral 
is presented in the Appendix. 

A typical airfoil surface is represented by a number of 
chordwise panels (N) over each of which the doublet strength 
(K n ) is assumed to be a constant. Kt can be expressed as a 
function of using Eqns. (3), (5) and (6) s 
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+ i^P N H- Jc (2 - 
e 11 '^N + 2iwD M 
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where 0^ is the semi-width of the last panel, N. For chosen N 
values ‘ of X at the collocation points of each of the N strips, 
using Eqns. (7) and (8), a set of N linear equations are ob- 
tained. For the known boundary conditions, W(X n ), the solution 
of these linear equations yields the K distribution from which 
the airloads can be computed. 

The above analysis is valid for any arbitrary motion where 
W and K value.- *re time-dependent. However, for aeroelastic 
analysis applications, it is a standard practice to assume a 
simple harmonic motion. For a two degree-of-f reedom simple 
harmonic motion in the plunge and pitch of the airfoil. 
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where Z* (** z / i , assumed to be positive downward) and a* are 
the amplitudes of the airfoil motion due to blade bending and 
torsion, respectively. The boundary condition at any colloca- 
tion point, i, then becomes 


W ( X , T) 


£iu>z' ♦ (1 ♦ iu»x^) a' J e luJ,r 
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where 


WfX^ = iwz' + (1 ♦ iwXj^) a'. 

Similarly, for a simple harmonic motion, K * where 

K' (* K ' ( X) ) is a function of x only. From these conditions, 
Eqns. (7), (8), (9) and (10) can be replaced with quantities 

that are fuir.'ons of X only, and the doublet distribution 
(amplitudes i be obtained for the assumed boundary condi- 
tions give oy Eqn. (11). 

The unsteady lift and moment coefficients are then expres- 
sed as 
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The parentheses, (•)* in Eqns. (14) and fl5) can be replaced 
either by z or a and the appropriate boundary conditions are 
applied . 


TWO-DIMENSIONAL FLUTTER ANAL V S IS 


Flutter is an oscillating instability produced by the nega- 
tive damping terms of the aerodynamic forces at speeds known as 
critical flutter speeds. The classical type of flutter is 
associated with potential flow and it is almost universal to 
maxe the assumption of linearli ty , according to which the forces 
brought into play by any small deviation and of its time deriva- 
tive The flutter problem can be reduced to two basic problems-- 
mechanical and aerodynamic. The first involves consideration of 
the motion of a flight structure as a continuous vibrating 
system acted on by external forces and internal damping, which 
then becomes one of writing the equations of motion for such 
a system. The second basic problem is that of determining the 
nature of the oscillatory aerodynamic forces, which are indepen- 
dent of the static forces that are needed to maintain the system 
in an equilibrium position. 

Equations of Motion 

Although an actual flight vehicle shouTd, strictly speaking, 
be considered a single elastic unit, it becomes necessary from 
the engineering viewpoint to make certain simplifying assumptions. 
Theodorsen (Ref. 18) and Theodorsen and Garrick (Ref. 19) ap- 
proached the flutter problem by considering a wing of infinite 
aspect ratio, moving with small oscillatory amplitudes at con- 
stant velocity. Thus the aerodynamic forces can be determined 
by considering the problem as one of two-dimensional flow. 
Furthermore, instead of the actual distributed mass and geometric 
properties of the wing, Theodorsen and Garrick assumed at that 
time that the results for the actual wing could conservatively 
be obtained by considering the motion of unit span of the wing 
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at some representative position; this has been often chosen 
at the 3/4 semispan. This approach is termed the two-dimen- 
sional flutter problem. 

To simplify the problem further, the issumption is made that 
the actual motion of a rotor blade system can be considered a 
combination of fundamental bending and fundamental torsion about 
the elastic axis. The system can then be replaced by an equiva- 
lent system containing an airfoil section of unit span restrained 
by springs against independent vertical (bending) motion and 
torsion. Then the equations of motion for a two degrec-of- 
freedom (bending torsion) system of an airfoil are 

M*2 + Sci + K z lZ - -pU 2 t <C tz Z + C £ a) (16) 


+ la + K a « pU 2 t 2 (C Z + C a) (17) 

a mz ma 


where , 


M 


l 

Z 

a 


K. 


K 


a 



dm = blade mass per unit span, 


r dm = blade static moment about elastic axis 
per unit span, 

r : dm = blade mass moment of intertia about 
elastic axis per unit span. 


= semi-chord length of the blade. 


= flapping (bending) motion, 

= pitchinj motion, 

= generalized bending stiffness of the blade, and 


= generalized torsional stiffness of the blade. 
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Assuming simple harmonic motion, 


z - z'e^ - z'e i “ T 

a » a'e^* - a'e i “ T 

and separating the aerodynamic derivatives 
ary parts, 


(18a) 

(18b) 

into real and imagin- 
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C = C + iC 
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(19a) 

(19b) 

(19c) 

(19d) 


Eqns. (16) and (17) reduce to 

(AK z - M£p 2 + pU 2 iC- z )Z* + (pU 2 ZC ia - Sp 2 ) a' = 0 

- ( pU 2 £ 2 C mz + S fcp 2 ) Z 1 + (K a - pU 2 p 2 C ma - Ip 2 ) a ' = 0. 


( 20 ) 

Equations (20) are a set of two linear, homogeneous equa- 
tions in Z ' and a' . For this system to have a nontrivial solu- 
tion, the coefficient determinant of Z' and o' must be equal to 
zero . 
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This determinant is referred to as the flutter determinant and 
the solution gives the flutter frequency. However, a direct 
solution cannot be obtained since the aerodynamic derivatives 
are functions of the reduced frequency. Since the aerodynamic 
derivatives are complex quantities, the determinant can be 
expressed in two parts, real and imaginary. After expandi r j 
the flutter determinant and substituting appropriate values 
of Eqns. (19) , the real and imaginary parts of the flutter 
determinant become, 
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where p = 1/p 2 . If the real and imaginary parts are set to zero, 
solutions p, and p 2 can be obtained for the real part, and p 3 
for the imaginary part. For an assumed value of w, if one of 
the solutions, p 3 or p ? , of the real part is equal to the solu- 
tion, pj , of the imaginary part, then this w corresponds to 
the flutter frequency. However, the aerodynamic derivatives 
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are functions of the Mach number and the reduced frequency. 

The flutter problem can only be solved after the aerodynamic 
derivatives are evaluated. It is necessary to assume a Mach 
number and a reduced frequency and test whether flutter occurs 
at these values. If the test results are negative, then it is 
necessary to iterate on reduced frequency until a flutter case 
occurs at these values. 


DESCRIPTION OF THE COMPUTER PROGRAM 


The main program incorporates the iterative scheme for 
a two-dimensional bend ing-torsion flutter analysis as described 
in the previous section. The geometric and structural data of 
a reference section of a wing of a rotor blade are read in as 
input information to the program. For any assumed value of 
reduced frequency (to) , the unsteady airload and moment coeffi- 
cients are obtained by calling one of the two subroutines, 
either AEROCC or AEROTH. Then in the main program the flutter 
determinant is formulated, the solution of which yields two 
real roots and one imaginary root. If one of the two real 
roots is not equal to the imaginary root, then the assumed 
reduced frequency does not correspond to a flutter condition 
and another value of w is assumed, and the process is repeated 
iteratively until the flutter frequency is obtained. From the 
flutter frequency, the flutter velocity and the Mach number are 
computed . 

Subroutine AEROCC is the unsteady airload prediction pro- 
gram for circulation-control airfoils. The numerical lifting 
surface method described in a previous section was adopted. 

The boundary conditions are assumed to correspond to a simple 
harmonic motion of the airfoil in plunge and pitch, with the 
fixed jet blowing momentum coefficient as a parameter. For a 
specified Cj, the steady airload and moment coefficients are 
computed as a function of reduced frequency and are used in the 
main program in evaluating various terms in the flutter equations. 

Subroutine AEROTH is the unsteady airload prediction 
program based on Theodorsen's hard flap solution for incompres- 
sible two-dimensional flows (Ref. 18). The theory and the 
computational procedure adopted is not presented in this report 
since it is available and well presented in the literature. 

The flap/chord ratio and the reduced frequency are the para- 
meters used in the program The program involves the genera- 
tion of Hankel functions as a function of to. The purpose of 
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this subroutine is to provide a cross check and also an alter- 
nate means to compute airloads in place of a circulation-control 
airload prediction. 

The input/output specifications and other relevant infor- 
mation for the computer program are included in the Program User's 
Manual . 


RESULTS AND DISCUSSION 


The numerical procedure adopted for the computation of 
airloads on circulation-control airfoils was based on a lifting- 
surface method in which a thin airfoil is divided into N strips 
over each of which the doublet strength is assumed to be a con- 
stant. The computed results for a typical case of an oscillating 
jet are shown in Table 1. From these results, it can be seen 
that the program yields a convergent solution even for the case 
in which the airfoil is represented by only 20 panels. 

Fiqures 2 and 3 present the lift response of an oscillating 
jet and its comparison with several theoretical and test results. 
The oscillatory jet deflection angle (amplitude) was assumed to 
b 6°, based on test results to correspond with Cj = 0.1 under 
steady flow conditions. Figure 2 presents the comparison of 
the lift amplitude ratio (|Cl|/|Cl 0 |) variation with the reduced 
frequency (w) . As can be seen from this figure, the amplitude 
ratio decreases as the frequency increases for most of the 
theoretical and test results, including the results of the 
present study. However, for the test results of Takeuchi (Ref. 

5) and the theoretical solution of Spence (Ref. 7) , the ampli- 
tude ratio increases as w increases. Takeuchi used false walls 
in his tests and the discrepancy may be due to these walls and 

tunnel interference. Spence's solution was based on a weak jet 

approximation and is valid for only very weak jet momentum 
coefficients. The present results are in good agreement with 
those of Theodorsen's hard flap solution (Ref. 18), and are in 
reasonable agreement with Potter's theory (Ref. 8) and the test 
results of Simmons and Platzer (Ref. 3) and Kretz (Ref. 6). 

Figure 3 presents the comparison of the phase lag angle 

with a). Once again Takeuchi 's and Spence's solutions are com- 
pletely out of phase from the other results. The present re- 
sults compare reasonably well with Theodorsen's and Potter's 
solutions for w values less than 1. They deviate substantially 
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Table 1. Lift Amplitude Ratios and Phase Angles for an 
Oscillating Jet (Cj ** 0.14 and Jet Deflection 
Angle = 7.3°). 

NUMBER OF PANELS 


20 24 30 


C _L_ 

# 

C L 


C L 

C Lo 


C Lo 


C Lo 


0.05 

0.923 

-0.126 

0.923 

-0.126 

0.923 

-0.124 

0.10 

0.857 

-0.176 

0.857 

-0.175 

0.858 

-0.173 

0.15 

0.804 

-0.197 

0.805 

-0.196 

0.806 

-0.194 

0.20 

0.762 

-0.205 

0.763 

-0.204 

0.764 

-0.201 

0.25 

0.728 

-0.204 

0.728 

-0.203 

0.730 

-0.200 

0.30 

0.700 

-0.193 

0.701 

-0.197 

0.702 

-0.193 

0.35 

0.677 

-0.190 

0.678 

-0.188 

0.679 

-0.1:4 

0.40 

0.658 

-0.179 

0.658 

-0.177 

0.659 

-0.173 

0.45 

0.641 

-0.167 

0.642 

-0.165 

0.643 

-0.161 

0.50 

0.627 

-0.155 

0.627 

-0.153 

0.628 

-0.150 

0.55 

0.615 

-0.143 

0.615 

-0.142 

0.616 

-0.138 

0.60 

0.604 

-0.131 

0.604 

-0.130 

0.605 

-0.126 

0.65 

0.594 

-0.120 

0.595 

-0.119 

0.595 

-0.115 

0.70 

0.536 

-0.110 

0.586 

-0.109 

0.587 

-0.105 

0.75 

0.578 

-0.099 

0.579 

-0.099 

0.579 

-0.096 

0.80 

0.572 

-0.090 

0.572 

-0.090 

0.572 

-0.087 

0.85 

0.566 

-0.082 

0.566 

-0.082 

0.566 

-0.080 

0.90 

0.561 

-0.074 

0.561 

-0.075 

0.561 

-0.073 

0.95 

0.557 

-0.067 

0.557 

-0.068 

0.557 

-0.066 

1.00 

0.553 

-0.060 

0.553 

-0.062 

0.553 

-0.061 

1.05 

0.550 

-0.054 

0.549 

-0.057 

0.549 

-0.057 

1.10 

0.547 

-0.049 

0.547 

-0.052 

0.547 

-0.053 

1.15 

0.545 

-0.045 

0.545 

-0.048 

0.545 

-0.049 

1.20 

0.543 

-0.040 

0.543 

-0.045 

0.543 

-0.047 

1.25 

0.542 

-0.037 

0.542 

-0.041 

0.542 

-0.044 

U) “ 

Reduced 

Frequency 

(pi/U) 

<J> - Phase 

Angle (radians) 


C T - Lift Coefficient 

L 

C T - Steady Lift Coefficient (w = 0) 
Lo 


15 


Lift Amplitude Ratio ( 


2. Or 


‘Spence (C J =0 . 1) 


u 


Present Analysis 
Other Theories 


1.0 


0.5 


Potter (C =0. 1) 

J 



Theodorsen's Hard 
Flap Solution 


Experiments 

A - Takeuchi (C =0 . 2 ) 

J 

□ - Kretz (C=0. 2) 

O ~ Simmons and Platzer 
(Cj=0 . 14 ) 


0.5 1.0 1.5 

Reduced Frequency (w) 


Figure 2. Comparison of Results for Lift Amplitude 
Ratio of an Airfoil with an Oscillating 
Jet at the Trailing Edge. 
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Figure 3. Comparison of Results for Lift Phase Angle 
of an Airfoil with an Oscillating Jet at 
the Trailing Edge. 
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front tho rest reaulfs as well oh from tho theoretical solutions 
for higher values of <• . However, note that for a typical classi- 
cal flutter analysis, only lower values of reduced frequency 
(u> < 0.5) are of interest. Also, the test conditions of the 
several experimental results are different and, hence, a good 
correlation is not really expected. 

The unsteady airloads program is computationally very 
efficient and it takes only about 0.02 cp seconds on a CDC 7600 
for one value of w. This is a very desirable feature, since 
a typical flutter analysis involves an iterative scheme re- 
quiring the computation of unsteady airloads for several values 
of 0) . 


The Flutter Program is checked out by applying it to a CCR 
blade case. The following data at 0.75 radius section of a CCR 
blade are taken from Reference 20: 


Reference chord length 

Mass per unit span 

Elastic axis location 

Moment of inertia per unit span 

First static moment per unit span 

First bending frequency 

First Torsional frequency 

Alternate torsional frequency 


1.467 ft. 

0.3737 slugs/ft. 

0.35c 

0.0776 slugs/f t . 2 /f t . 
0.0547 slugs/f t. 2 /ft . 
4.8 cps 
44.5 cps 
15.0 cps 


The flutter analysis was performed for this airfoil assuming 
Cj = 0.1 at a jet (or flap) deflection angle of 6°. The computed 
flutter speeds using each of the two airload programs, AEROCC and 
AEROTH, are tabulated below. 


Altitude 

Torsional 

Flutter Speed 

(ft. /sec. 


Frequency 

(cps) 

AEROCC 

AEROTH 

Sea level 

44.5 

1,248 

1,252 

10,000 ft. 

44.5 

1,448 

1,454 

Sea Level 

15.0 

376 

379 

10,000 ft. 

15.0 

435 

4 39 


This example once again demonstrates that the circulation- 
control airloads program developed under the present study is in 
agreement with Theodorsen's hard-flap solution. Also note that 
the airloads computed in this program were based on the incompres 
sible flow assumption, and hence the higher flutter speeds 
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obtained using this program are in gross error and are beyond 
the range of the Program's validity. 


CONCLUSIONS 


A fast and efficient computer proqram for predicting the 
unsteady airloads in incompressible flow was developed and 
correlated. The results for CC airfoils under harmonic motion 
are in good agreement with those of Theodorsen's hard-flap 
solution. A two-dimensional bending-torsion flutter computer 
program was developed and used by performing a flutter analysis 
of a typical CC rotor blade section. 

In the present analysis, the compressioility and the rotor 
wake effects are not taken into account. As the flutter Mach 
number increases, it is important to include the compressibility 
effects, and usually they tend to decrease the flutter speed. 

In a typical hovering rotor, the wakes of the other blades are 
directly underneath the reference blade and may have a signifi- 
cant effect on the flutter speed. The rotor wake effects can 
be incorporated in the present numerical lifting-surface method 
by adopting the mathematical model of Loewy (Ref. 21). The 
compressibility effects can be incorporated by adopting the 
analysis of Jones and Rao (Ref. 22) . 
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APPENDIX: EXPANSION OF WAKE INTEGRAL 
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